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General properties of global covariance matrices representing bipartite Gaussian states can be 
decomposed into properties of local covariance matrices and their Schur complements. We demon- 
strate that given a bipartite Gaussian state pi2 described by a 4 x 4 covariance matrix V, the Schur 
complement of a local covariance submatrix Vi of it can be interpreted as a new covariance matrix 
representing a Gaussian operator of party 1 conditioned to local parity measurements on party 2. 
The connection with a partial parity measurement over a bipartite quantum state and the determi- 
nation of the reduced Wigner function is given and an operational process of parity measurement is 
developed. Generalization of this procedure to a n-partite Gaussian state is given and it is demon- 
strated that the n— 1 system state conditioned to a partial parity projection is given by a covariance 
matrix such as its 2 x 2 block elements are Schur complements of special local matrices. 

PACS numbers: 



INTRODUCTION 



Several aspects of the Quantum Information theory rely on the ability to prepare and manipulate the quantum state 
of a given physical system. Obviously this ability is related to the available physical operations on the actual state 
encoding. For quantum protocols whose state encoding is spanned by a two dimensional Hilbert space (qubits) one 
has to be able to completely describe the actual physical operations acting directly on the physical system codifying 
state space. Thus any physical operation (including measurements) is defined through completely positive maps. For 
continuous variable encoding as given by Gaussian states [33| physical operations given by positive maps can or cannot 
keep the Gaussian character of the system state. Those operations that map every Gaussian input state into a Gaussian 
output state are called Gaussian operations [ij . For unitary operations acting on a single system Gaussian state, for 
example, the positive map is completely described by a local Sp(2, R) symplectic (canonical) transformation over the 
system's covariance matrix. Conversely, as a consequence of the Stone-von Neumann theorem, to every symplectic 
transformation over the covariance matrix there exists a unique unitary operation acting on the corresponding state 
space. Those operations correspond to the transformations induced by linear (active and passive) optical devices 
(beam-splitters, phase-shifters, and squeezers). It is thus not surprising that continuous variable quantum information 
protocols have been implemented with a Gaussian state encoding and Gaussian operations 

The whole set of unitary operations does not describe the most general transformation on quantum system states. 
A general evolution map of a quantum system including coupling to ancillas (or reservoir) and measurement is given 
by quantum operations [4J. In this formalism the map 

- . *W (1) 



Tr[e{p)\ 



connects the initial (input) to the final (output) state. The quantum operation e is a linear, trace decreasing map that 
preserves positivity. A quantum operation s satisfying complete positivity is written as e(p) = JT,j AjpAp where Aj 
is a set of system operators, which must satisfy J2j AjAj < I. Generally speaking, "interactions" with an ancilla (or 

reservoir) satisfy the trace-preserving condition Y). A^A-j = I, while quantum measurements satisfy Y]j AjAj < 1. 
Those operations which are trace preserving and irreversible are sometimes named as quantum channels. Gaussian 
quantum channels [B| have the additional feature to keep the Gaussian nature of the systemfs) quantum state. 

Several authors have described Gaussian quantum channels and operations throughly [E[l,[l,0,HS[l3]- Regardless 
the specific interest in each one of those references, one common interesting feature observed is that the Schur 



complement ll| of square matrices representing Gaussian states covariances embodies a manifestation of a physical 



operation when considering partial projections and trace operations onto Gaussian states 0, Q ■ However the process 
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that delivers the Schur complement of the covariance matrix itself for an input two-mode bipartite Gaussian state has 
not been discussed yet. As we show in this paper, this is only possible when non-positive operations are allowed. The 
purpose of the present work is to derive this process for an arbitrary input two-mode bipartite Gaussian state and 
generalize it for the case of an partite (n-modc) Gaussian state. Particularly we show that the Schur complement 
of one of the local covariance matrices of a bipartite Gaussian state appears as the covariance matrix describing a 
Gaussian operator, which embodies a parity measurement over the other mode of the input state. 

We begin by reviewing in Sec. 2 some properties of bipartite Gaussian states |l2J and we describe the decomposition 



of positivity for a covariance matrix through the Schur complement structure. In Sec. 3 we discuss examples of 
Gaussian quantum operations, consisting of the vacuum state projection on one mode of a two-mode bipartite input 
state, and the partial trace operation. In Sec. 4 we develop our central result, namely the covariance matrix that 
represents the Schur complement of one of the input subsystem covariance matrix. In Sec. 5 we describe the parity 
measurement process and the physical properties behind the mathematical quantity as given in Sec. 4. In Sec. 6 
we outline how such a Schur complement of a local covariance matrix can be employed for entanglement and global 
purity quantification. In Sec. 7 we provide a generalization of this local parity measurement for a Gaussian state with 
n-modes in terms of an output covariance matrix. Finally Sec. 8 concludes the paper. 



TWO-MODE BIPARTITE GAUSSIAN STATES 



Any two-mode bipartite quantum state, pi2> is Gaussian if we can write it as [34 

Pl2 = [dz e ^Ea e -|ztVz_ 



(2) 

Alternatively, if its symmetric characteristic function is given by \( z ) — Tr[D(z)p] — e _ i zt ^ z , where D{z) = e z ^ a 
is a displacement operator with zt = (z*, z\, z 2 , z±), at = (aj, ai, a\, a±), being ai(aj) and a2(c4) the annihilation 
(creation) operators for modes 1 and 2, respectively. Remark that here, as in many integrals through this paper the 
following notation is employed: dz = dz\dz 2 , where dzi = d 2 Zi = dIm{zi)dRe{zi) for i = 1,2. Also, 
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and V is the Hermitian 4x4 covariance matrix describing all the second order moments Vij = (— 1) I+J ({v i: vj})/2, 
where v 2 , v 3 , u 4 ) = (ai, a\ , a 2 , a\), given by 
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where Vj and V 2 are 2x2 Hermitian matrices containing only local elements while C is the correlation between the 
two parties. Furthermore any covariance matrix must be positive semidefinite (V > 0) and for the operator (T5|) to 
represent a physical state the inequality 



V+iE>0, 



(5) 



must hold, which is nothing but the fundamental uncertainty principle. 

Those general positivity criteria can all be decomposed into properties of the block matrices Vi,V 2 and C. A 
convenient decomposition procedure is given by the Schur complement for non-singular block matrices 35| . In this 
way for the covariance matrix ((4]) , V > if and only if 



V, > 



(6) 



and the Schur complement of V 2 , given by 



S(V 2 ) = Vi - CVg X C f > 0. 



(7) 
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Similar properties have to be satisfied for the generalized uncertainty principle 



and 



v 2 + iz > 0, (8) 



V! + -ZJ - C (v 2 + -Zj Ct > 0. (9) 

The emergence of the Schur complement structure and its importance can be fully appreciated through the well known 
determinant decomposition property. Given V as described above, 

detV = detVidetCVa-C+V^C) 

= detV 2 det(Vi - CV^C*). (10) 



Thus, it is immediate that any operation that yields Vi and S(Vi) (or V2 and S(Va) by symmetry) can be employed, 
in priciple, to quantify some physical properties of the global covariance matrix V. This is the case for entanglement 
and global purity of a bipartite system, as discussed in details in Ref. [HI], whose central results are commented latter. 



GAUSSIAN OPERATIONS 

Let us remark on two examples of Gaussian operations, which were previously considered in details in Ref. [H, 0, 0] 
to show that Schur complements of block matrices representing Gaussian states covariances embodies a manifestation 
of a physical operation when considering partial projections onto Gaussian states. 

The first example considered is a vacuum state projection over the mode 2 of a two- mode Gaussian state p\ 2 : 

<r[ 0) =7> 2 {|0) 2 (0| Pl2 }. (11) 
Writing the Gaussian state form of the one mode vacuum state according to equation Q in pip , we have 



a[ 0) = Tr 2 {dz 2 e-i z ^e- z ^ Pl2 }. (12) 



Now replacing the Gaussian form for p\ 2 ([2]) and evaluating the trace, we have that 



det (V 



1 f dzi ^tZa^zt^-Cf^irCtjz! 

V 2 + il) J 



(13) 



with dzi = d 2 z\ = dIrn(zi)dRe(zi). We notice that the resulting state has also a Gaussian form with covariance 
matrix given by (|14[) 

lf } = Vi - C (V 2 + llj Ct. (14) 

Here I is a 2 x 2 identity matrix. The covariance matrix of the reduced state, (fT4"f is the Schur complement of the 
matrix 

m ' = (c! Vj + ii) < 15 > 

in relation to V 2 + |l [l|. 

In the second example we consider a simple partial trace operation: 

Pi = Tr 2 {pi 2 }. (16) 
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In this case it is obvious that p\ is still a Gaussian state and that it will be given, in terms of the representation [2j as 

Pl = J dz ie z l Z ^e-i 2 l Y ^. (17) 

Now, the covariance matrix of the reduced state is given by the Schur complement of matrix V2 as given by the 
matrix 

M *=(^v 2 )- ( 18 ) 

EMERGENCE OF THE SCHUR COMPLEMENT STRUCTURE OF LOCAL COVARIANCE MATRICES 

As we noted in the previous section, projection operations followed by partial trace operations are manifested as 
Schur complements of 4 x 4 matrices. However none of those matrices is the covariance matrix of the input bipartite 
Gaussian state since the vacuum state covariance matrix must be inserted, or the correlation matrix C neglected. 
Therefore, we need a physical measurement process that is able to yield the exact Schur complement of one of the 
local covariance matrices. In this section we describe the mathematical operation which is able to realize this task 
and in the next section we describe the physical operation behind it. 

Proposition 1: Given a bipartite Gaussian state P12, the covariance matrix Ti describing the Gaussian operator o\ 
of mode 1 conditioned to a parity projection over the mode 2, 

ox =7> 2 {e^ Q2 pi 2 }, (19) 

is given by the Schur complement of the covariance matrix V of the input bipartite state p\2 in relation to V 2 : 

Ti = Vi - CV^C*. (20) 



Proof: We can write the parity operator as an integral over the displacement operator [22], |24j], in this case the 
equation (TlT)]) can be rewritten as 

vi = ^Tr 2 ^J cfe 2e -^ Za V 12 j, (21) 

with ds2 = d 2 S2 = dIm{s2)dRe(s2) ■ As pi 2 represents a Gaussian state we can write it in the form ([2]) and replacing 
this Gaussian form in (f2"Tj) we have 

= l Tr2 || ds2e -st Z a 2 | rfzi | rfz2e ztZa le ztZa 2e -iztVz| 

= \jds 2 J dz,J dz2Tr 2 {e^ s ^ Z ^}e-i z l Zz - 

xe ztZa le -iztVz e iztZs 2i (22) 



but 



therefore 



Tr 2 j e (z 2 -s 2 )+Za 2 j = e -||* 3 - S2 |^(2)( Z2 _ ^ (23) 



izi'Vz 



o-i = ~ J dz ie z i Za > J dz 2 e~3' 

J dz ie < Z ^e-^- C ^ 1C >\ (24) 



1 1 



2 VdctV 2 

where we notice that the resulting state has also a Gaussian form with covariance matrix given by (|20 
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EXPECTED PARITY AND PROBABILITY DISTRIBUTION FUNCTIONS 

In the previous section we have given a mathematical connection of the Schur complement of a local covariance 
matrix V 2 with the parity operation e l7ra 2 a2 over mode 2 of a two- mode Gaussian state and further reduction to 
mode 1. To completely understand the physical meaning of such an operation we must develop a theory for parity 
measurement detection. 

Indeed, parity operation properties have received considerable attention recently, in connection to experimentally 
accessible measures, for description of quasiprobability distribution functions for single and entangled systems and 
for quantum information and computation proposals wi EE E3, EE EE H3] ■ Moreover the average parity of a given 



quantum state is related with its Wigner function |2l| at the origin of the phase space and nowadays it can be 
experimentally determined for radiation fields through photocounting experiments [23 . [23} , or in microwave cavity 
quantum electrodynamics experiments [U l25| 

Bellow we give a description of the parity measurement process. Given a prior quantum state p, the post-selected 
state conditioned to a parity (p = 1, — 1) measurement is given by 

»-T%fry (25) 

with probability P p = Tr {J\f p p}, where the super-operator J\f p , p = 1, — 1 is given by§S6j 

M p -=J2\n)(n\-\n)(n\, (26) 

nGp 

i.e., it indicates projections over even or odd Fock states \n). For p = 1 the sum runs over even natural numbers, while 
if p = — 1 it runs over odd natural numbers. The probability for the occurrence of the two events write independently 
as 

P 1 =Tr{Af 1 p}=Y^ ( n \p\n), (27) 
P_i = rr{JVlip} = X>Mn> s (28) 

n a dd 

and it can be evidenced that P\ + P_i = 1, as it should be. The average parity, p, is then simply given by 

P = J2P P P=T, H/>l»> ^ E H^> = Tr{e^ a p) . (29) 

P Tl dd 

The average parity was recognized independently by Grossmann [27j ] and by Royer [28| to be proportional to the 
Wigner distribution function at the origin of the phase space: 

W(0) = 2 p, (30) 

The values of this function in other points of the phase-space can be achieved by performing displacements over the 
input state 0, 0, HI , such that 



W(a) =2j2(-l) n {n\D( a )prf(a)\n), (31) 

n=Q 

where D(a) is the displacement operator. 

This approach can be extended for bipartite systems as given by the density operator pi2- The joint post-selected 
state P12 conditioned to a parity (p — 1,-1) measurement over the mode 2 is given by 

$ - rf" 12 y (32) 

Tr{N p P12} 

with probability P p — Tr 12 {Af p P12}, and the operation Af p , now reads as 

7V p - = Ii<g)^|n)<n| 2 »<«| 2 . (33) 



6 



Now the probability for the occurrence of the two events write as 

Pi = Tr 2 {Wi p 2 } = H<°2k>2, (34) 

P-x = Tr 2 {M-i p 2 } = ^2{n\ P 2\n) 2 , (35) 

"odd 

where p 2 — TVi{pi2}. The average parity of mode 2 is related to its Wigner function as 

W 2 (0) = 2 p = 2 Tr 12 {e*™!^ p 12 j , (36) 

In view of this, the significance of the operator <j\ is clearly identified as the difference between the states of mode 1 
conditioned to projections on the even Fock subspace and on the odd Fock subspace of mode 2: 

°i = ^2 (n\pi2\n) 2 -^2(n\p 12 \n) 2 , (37) 

"raen "odd 

i.e., the difference between the states of mode I conditioned to even and odd parity measurements on mode 2 out of 
an ensemble of identicaly prepared pi 2 states. 

Finally, an interesting relation follows from Eq. (). It is immediate to check that the Wigner function of mode 2 at 
the origin of the phase-space can be inferred by 

W 2 (0) = 2Tr 1 {a 1 }. (38) 

In this case, this result becomes very interesting, for it means that in a bipartite Gaussian state, the achievement of 
the Wigner function at the origin of the phase-space of subsystem 2, is the trace of the Gaussian operator of subsystem 
I (7i, as above, and vice versa. The resulting covariance matrix representing <t\, will be given by (|20| . 

Also, from Eq. ([36]), we can derive a relationship between the Wigner function at the origin of phase space IF2(0) 
of one mode Gaussian state with covariance matrix V2 and the Schur complement S(n 2 + i) of one diagonal clement 
of V2, which in this case is a scalar, being 

W 2 (0) = - 1 (39) 

^{n 2 + \)S{n 2 + ±) 



where 



F , 1 ITO2I 2 



S(na + r)=n 2 + --^ T . (40) 
2 2 n 2 + ^ 

Moreover, we note that the Wigner function of a mode at the origin of the phase-space, and thus S(n 2 + i) is 
directly related with the I 2 = (n 2 + i) 2 — \m 2 \ 2 element of the four invariant set of the Sp(2, R) <g) Sp(2, R) group 0]: 
h = detVi, I 2 = detV 2 , h = detC and h = Tr[ViZCZV 2 ZC t Z]. 
Another way to write o\ is by expanding it in a coherent states basis, 

Px% = yd 2 /3id 2 /3 2 P(/3i,/32)|/3i,/32)(/?i,/32| (41) 

where P{j3\,[3 2 ) is the Glauber P function of the joint system p\ 2 . Replacing (|4i"j) in ([19]) and using ([36]) we have 

a x = [ d 2 /?i|/3i)(/3i| / d 2 p 2 W c {Q)P{Pi,p2), (42) 



where Wc(0) = e 2 '' 32 ' , is the Wigner function of the coherent state \(3 2 ) (that spans the base of mode 2) at the origin 
of phase space. Thus by defining 

V(Px) = [ d 2 p 2 W c (0)P(Px,(3 2 ) (43) 
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we have 

tri= J dWiXftTOi). (44) 

When there is only a reduction of the system without association to a measurement in one subsystem, the P function 
of the reduced state can be obtained by integrating the bipartite state P function, P[j3\,j32), over the traced mode 
variable. But in (|41[) , the P function is associated with a weight in the form of a Gaussian function, where it shows 
that the values of the P function of <j\ over the variables nearly to the origin of the phase space are more important. 
The influence in mode 1 by a parity measurement in mode 2 can also be noted if we compare the covariance matrix 
pop with the respective matrix Vi of the reduced operator p\ = Tr2{pi2}, indicating that the parity measurement 
insert global properties in the local terms of mode 1. 

ENTANGLEMENT CHARACTERIZATION OF A TWO-MODE BIPARTITE GAUSSIAN STATE 

As we saw in the previous sections, the achievement of the Schur complement structure of a two-mode Gaussian state 
covariance matrix is related with local parity measurements and the Wigner function at the origin of the phase-space. 
However, this result has a main physical interpretation related to the entanglement characterization of a two-mode 
Gaussian state via a protocol based on local operations and classical communication (LOCC) as discussed in Ref. 
0. 

The protocol consists basically in the attainment of all local sympletic invariants (Ji, Jg, I-Z3I and J4), discussed in 
the previous section, via only LOCC and through Proposition 1 of this paper. 

Suppose that Alice and Bob share many copies of a two-mode Gaussian state, where the subsystens related to 
mode 1 are given to Alice and those related to mode 2 to Bob. The quantities I\ and I2 can be locally determined 
by the reconstruction of the submatrices Vi and V2 or by purity measurements of the mode 1 and 2, performed by 
Alice and Bob, respectively. If Alice is the one who want to achieve the entanglement characterization, Bob has to 
send his matrix or purity outcome to Alice via a classical communication channel. I-Z3I and I4 are determined using 
the Proposition 1, where essentially Alice has to construct the covariance matrix Ti (|20|) from the diference between 
the correlation matrices of the mode 1 subsystem conditioned to projections onto odd and even states of the mode 2 
subsystem. For that, Bob has to perform a parity measurement on each subset component letting Alice to know to 
which copy does that operation corresponds and the respective outcome, i.e., even parity (eigenvalue 1) or odd parity 
(eigenvalue -1). Alice then separates her copies in two groups, the even (e) and the odd (o) ones, each one containing 
the copies conditioned by an even and odd parity measurement on Bob's copies respectively, allowing her to obtain 
the correlation matrices Vi e and Vi D related to each group. Consequently, the matrix Ti is obtainned by subtracting 
the odd correlation matrix from the even one Ti = Vi e — Vi Q . With this three matrices (Vi, V2 and Ti) in hand, 
Alice is able to completely characterize the Gaussian state's entanglement content as well as its purity without any 
global or nonlocal measurements. Two invariants are obtained by I\ = detVi and I2 = detV2, while the third one 
I-Z3I is calculated from the expression 



|/ 3 1 = ^2dct(V!-ri), (45) 



and the fourth quantity I4 by 



J4 = Ii/3+/3-/adetr 1 . (46) 

This four quantities allows, for example, decide whether or not the two-mode Gaussian states is entangled. For 
that, the Simon separability [29( is the relation used to perform this test, i.e., the state is not entangled if, and only 
if, 

hh +{\~ l J 3| 2 ) ~ J 4 > (h + J 2 )/4. (47) 
Moreover, for a symmetric state (Ii — I2), the entanglement can be quantified via the entanglement of formation 



(E f ) m\m- 
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Efipu) = / (z\fli + \h\-Vh + 2/i|/ 3 |) , (48) 

where /(x) = c + (x) log 2 [c + (x)] — c_(x) log 2 [c_(x)] an d c+ (x) = (x 1 / 2 ± x 1 / 2 ) 2 /4. For arbitrary two-mode Gaussian 
states (Ii ^ J 2 ) we can work with lower bound for Ef |30| or calculate its negativity or logarithmic negativity (32)]. 

Note that all the procedures adopted to achieve the four invariants were made without any type of global mea- 
surements or operations, allowing to characterize the two-mode Gaussian state's entanglement content via only local 
measurements and a classical communication channel. 



PARITY MEASUREMENT IN A iV-MODE GAUSSIAN STATE 

Now we can generalize the result of theorem 1 for the case of a ro-mode Gaussian state. 

Proposition 2: Parity measurement in the m mode of a Gaussian state of n-modes with covariance matrix 2n x 2n 
given by 

C 13 
V33 



/V u C 12 

c f 12 v 22 



V, 



4 



'13 



'23 



(n— X)n 



'(n— l)n 

v 

v nn 



(49) 



associate the state of n — 1 resulting modes, in such a manner that the resulting covariance matrix 2(n — 1) x 2(n — 1) 
is formed by 2 x 2 block matrices, located in line i and column j, given by 

(50) 

such that for i = j, Mjj = Vj is the covariance matrix of the reduced operator for the subsystem i, and for i ^ j , 
My = Cij are the matrices that represent the correlations between the n modes of the global system, noting that in 
this case = Mj • . 



T tj = M y - M^M-^Mt 



Proof: This proof is made by induction. As we have already derived how a parity measurement in one mode of 
a bipartite Gaussian state affect the other mode in terms of the covariance matrix (|20| . it is possible to derive a 
simmilar relation for states with 3 and 4 modes and we can verify that there is a standard form between the influence 
of a parity measurement with the resulting reduced covariance matrices, allowing to make a generalization in the case 
of a n-mode state. For a tripartite Gaussian state, P123, with covariance matrix given by 

Vi C12 C13 



V123 = I C T 12 V 2 C 23 I , (51) 

"l3 C L V 3 

the resulting bipartite Gaussian state will be conditioned to a parity measurement in the mode 3 of the global system, 



<7i2 = 2>{e ITra 3 a3 pi23}, such that the respective covariance matrix is 



2 — 



Vi - C13V3 ^t, Cia - C 13 V3 X C 



'12 



C 2 3V3- 1 c" 3 V 2 



3 ^23 

C23V3 c 23 



(52) 



Note that each block element of the matrix above are a Schur decomposition of an another matrix 4x4 and can be 
obtained by the relation (|50|) . 

With an analogous calculus for the case of a 4-mode state, the reduced covariance matrix conditioned to a parity 
measurement in mode 4 is given by 

Vi — C14V4 X C 14 C12 — Ci4V 4 1 C 24 C13 — Ci4V 4 1 C 34 

ri M =i cL-c^vr'cL v 2 -CmV7 1 cS„ c M - Cwvr^t I. (53) 
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ChaVa x C\a C 23 — CzaVa 1 Ct A V 3 — ChaVa 1 Ci 



In the same way of the case of a measurement over the states with 2 and 3 modes, we noted that the block elements 
of the matrix given by (f5"3")) can also be described by (|50p . 
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CONCLUSION 



In this work, we have investigated which operation over a two-mode bipartite Gaussian state delivers the Schur 
complement form of an input local covariance matrix. We discovered that parity measurements in one mode of a 
global bipartite two-mode state influences the other mode in such a manner that there exists a Gaussian operator 
whose covariance matrix has the form of a Schur complement of one of the local covariance matrices of the input 
state. This operator is given as the difference between the reduced state of one subsystem conditioned to even and 
odd projections on the other subsystem. As the parity measurement and the Wigner function are strictly related, it 
is possible to associate the Schur complement of a local covariance matrix with the process to achieve the Wigner 
function of one mode and, due to the invariance of the Gaussian properties by displacements, with the achievement 
of only one point of this function. At the origin, the Wigner function of an one mode Gaussian state is related with 
one element of the four invariant set of the Sp(2, R) <g> Sp(2, R) group and so is the Schur complement of a local 
covariance matrix. Moreover we have generalized the approach for a n-partite Gaussian state verifying that after a 
parity measurement in one mode, the n—1 system state has a covariance matrix with 2x2 block elements in a form 
of Schur complements of special block matrices. We believe that our findings have both conceptual and practical 
implications for the development of continuous variable protocols with Gaussian states [Hj]. 

We are pleased to thank G. Rigolin for valuable discussions. This work is partially supported by FAPESP and by 
CNPq. 



References 



[1] Eisert J and Plenio M B 2003, Int. J. Quant. Inf. 1, 479. 

[2] Braustein S L and Kimble H J 1998, Phys. Rev. Lett. 80, 869. 

[3] Furusawa A et al. 1998, Science 282, 706. 

[4] Kraus K 1984, States, Effects, and Operations (Berlin:Springer-Verlag). 
[5] Eisert J Wolf M M 2005, preprint |arXiv:quant-ph/050515i~j 
[6] Giedke G and Cirac J I 2002, Phys. Rev. A 66, 032316. 
[7] Eisert J, Scheel S and Plenio M B 2002, Phys. Rev. Lett. 89, 137903. 
[8] Fiurasek J 2002, Phys. Rev. Lett. 89, 137904. 
[9] de Oliveira M C 2004, Phys. Rev. A 70, 034303. 
[10] de Oliveira M C 2005, Phys. Rev. A 72, 012317. 

[11] Horn R A and Johnson C R 1987, Matrix Analysis (Cambridge:Cambridge University Press). 

[12] Englert B G and Wodkiewicz K 2003, Int. J. Quant. Inf. 1, 153. 

[13] Haruna L F, de Oliveira M C and Rigolin G 2007, Phys. Rev. Lett. 98, 150501. 

[14] Glauber R J 1963, Phys. Rev. 131, 2766. 

[15] Ozorio de Almeida A M and Brodier O 2004, J. Phys. A: Math. Gen. 37, L249. 

[16] Ozorio de Almeida A M, Vallejos R O and Saraceno M 2005, J. Phys. A: Math. Gen. 38, 1473. 

[17] Ozorio de Almeida A M, Preprint |quant-ph/0612029| 

[18] de Oliveira M C and Milburn G J 2002, Phys. Rev. A 65, 032304. 

[19] de Oliveira M C and Munro W J 2000, Phys. Rev. A 61, 42309. 

[20] Dodonov V V 2007, Phys. Lett. A 364, 368. 

[21] Wigner E 1932, Phys. Rev. 40, 749. 

[22] Banaszek K and Wodkiewicz K 1996, Phys. Rev. Lett. 76, 4344. 

[23] Banaszek K, Radzewicz C, Wodkiewicz K and Krasiriski J S 1999, Phys. Rev. A 60, 674. 

[24] Lutterbach L G and Davidovich L 1997, Phys. Rev. Lett. 78, 25477. 

[25] Bertet P et al. 2002, Phys. Rev. Lett. 89, 200402. 

[26] Cahill K E and Glauber R J 1969, Phys. Rev. 177, 1882. 

[27] Grossmann A 1976, Commun. Math. Phys. 48, 191. 

[28] Royer A 1977, Phys. Rev. A 15, 449. 

[29] R. Simon, Phys. Rev. Lett. 84, 2726 (2000). 

[30] G. Rigolin and C. O. Escobar, Phys. Rev. A 69, 012307 (2004). 

[31] G. Giedke, M. M. Wolf, O. Kriiger, R. F. Werner and J. I. Cirac, Phys. Rev. Lett. 91, 107901 (2003). 
[32] G. Vidal and R. F. Werner, Phys. Rev. A 65, 032314 (2002). 

[33] A Gaussian state is the one which can be fully described by its first and second moments, the latter of which are commonly 
collected in the covariance matrix. 



10 



[34] Throughout this paper we consider the first order moments set to zero with no consequence to the generality of the results, 
since this can always be performed with local displacements without altering the state entanglement properties. 

[35] Local covariance matrices, Vi and V2 are always non-singular due to the presence of the vacuum fluctuation term 1/2 in 
all their diagonal elements. 

[36] The construction of this quantum operation is derived inversely by recalling that the Wigner function at the origin of the 

phase space is the average of the parity operator. Indeed from that one would obtain M p - — (^2 nep \n) {n\j ■ (^2 nep \n) {n\j . 

The completeness relation for the Fock states followed by the fact that Af p - for all pratical purposes is always followed by 
the trace operation lead to the diagonal form of Eq. () convenient for the operator sum representation. 



